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Abstract 

X-ray line profiles represent a new way of studying 
the winds of massive stars. In particular, they en- 
able us to probe in detail the wind-wind collision in 
colliding wind binaries, providing new insights into 
the structure and dynamics of the X-ray-emitting re- 
gions. We present the key results of new analyses of 
high-resolution Chandra X-ray spectra of two impor- 
tant colliding wind systems, j 2 Velorum and WR140. 
The lines of 7 2 Vel are essentially unshifted from their 
rest wavelengths, which we suggest is evidence of a 
wide shock opening angle, indicative of sudden ra- 
diative braking. The widths of the lines of WR140 
are correlated with ionization potential, implying non- 
equilibrium ionization. The implications of these re- 
sults for the radio emission from these systems are dis- 
cussed, as are some of the future directions for X-ray 
line profile modelling of colliding wind binaries. 

1 Introduction 

Massive, early-type stars (M > 1OM ) have profound 
effects on their environments through their intense UV 
radiation fields, their powerful stellar winds (M ~ 
KT 7 -1(T 5 M Q yr" 1 , ~ 2000 km s" 1 ) and their 
violent deaths as supernovae. The winds from mas- 
sive stars provide ~30% of the mechanica l energ y in- 
put into the interstellar medium dAbbottL 119821) . and 
the mass-l oss also substantially af fects the stars' own 
evolution ( Chiosi & Maeder, 1986b . Thus, understand- 
ing mass-loss from massive stars is essential for un- 
derstanding feedback to the ISM, stellar evolution and 



galactic evolution. 



Most massive stars resi de (or r esided) 
or multiple systems dzinnecker . 2003). 



in binaries 
In such 

systems, the winds of the massive stars collide 
highly supersonically (Mach number > 100), and 
they are compressed and shock-heated to ~10 7 K. 
The shock-heated gas produces co pious X-rays (e.g. 
Stevens. Blondin. & Pollockl 1992), though there are 
observational signatures throughout the electromag- 
netic spectrum, including non-thermal radio emission 
from relativistic elec tr ons a c celerated at the shocks 
(e.g. bichler & Usov . 1993 : Dougherty & Williams, 
200(1 

Previously, studies of X-rays from colliding wind bina- 
ries (CWBs) have tended to focus on their broad-band 
spectral properties, due to the poor spectral resolution 
of the satellites then in operation (e.g. E/ AE ~ 20 for 
the ASCA SIS). Nevertheless, by comparing the vari- 
able ASCA spectrum of j 2 Velorum with hydrody nam- 



ical models of the system, Stevens et al. ( 1996) were 



able place constraints on the stars' mass-loss rates and 
wind velocities. This technique has more recently been 
applied to the hig h-resolution Chandra grat ing spec- 
trum of rj Carinae dPittard & Corcoranl.l2002h . 

The unprecedented spectral resolution offered by the 
gratings on board Chandra and XMM-Newton enables 
us to resolve line shifts and widths down to a few hun- 
dred km s , providing a unique probe of the dynam- 
ics of the wind-wind colllision. Furthermore, the X- 
ray forbidden-intercombination-resonance (fir) triplets 
from He-like ions provide diagnostics of the density 
and temperature of the X-ray-emitting plasma, as well 
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Figure 1 : The co-added first order MEG spectrum of 7 2 Vei. The spectrum has been binned up to 0.02 A. The 
absorption at iong wavelengths is due to the unshocked stellar winds. 



as potentially p roviding diagnostics of the UV radi- 
ation field (e.g. Paerels & KahJ 2003 ). Using these 
tools we can obtain new insights into the location, ge- 
ometry, structure and dynamics of the wind-wind col- 
lision region. 



2 Velorum 



7 



We have carried out a new analysis of an archived 65- 
ks Chandra High-Energy Transmission Grating Spec- 
trometer (HETGS) observ ation of the well-s tudied 
WR+O binary 7 2 Velorum (Skinne r et all 1200 ll : Hen- 
ley et al., in prep.). Vel is a double-lined 
spectroscopic binary of spe ctral type WC8 + 07.5 
( de Marco & Schmuta. 1 1999b whose orbit is well de- 
termined, with a period of 78. 53 ± 0.01 days , e = 



1997) 



0.326 ± 0.01, uw b = 68° ± 4° (Schmutz et al 
and i = 63° ± 8° Jde Marco & SchmutzLll999l 

The first order Medium Energy Grating (MEG) spec- 
trum of 7 2 Vel is shown in Fig. ^ The spectrum is 
dominated by strong emission lines from S (A « 4- 
5 A), Si (A » 5-7 A) and Mg (A » 7-9 A), with 
weaker lines from Ne (A w 9-14 A) and Fe (e.g. 
A « 2 A and A « 10 A). 

We have measured the line centroid shifts and the line 
widths by fitting Gaussian line profiles (plus a constant 
continuum) to a narrow wavelength region around the 



particular line of interest. Although colliding wind bi- 
naries are e xpected to exhibit a wide range of line pro- 
file shapes (Henl ey. Stevens. & PittardL 2003 ). Gaus- 
sians give good fits to the emission lines in 7 2 Vel and 
thus provide a good way of quantifying the line shifts 
and widths. 

The measured line shifts and widths are plotted against 
wavelength in Fig. 13 The lines are generally unshifted 
from their lab wavelengths. The mean shift is —52 ± 
16 km s _1 and the mean FWHM is 1240 ±40 km s _1 . 
Neither the line shift nor width is correlated with wave- 
length or ionization potential. 

In order to understand our measured shifts and widths 
in terms of the geometry of the wind-wind collision, 
we use a simple model in which th e wind-wind inter- 
action region is a conical surface (Liihr 1 Il997h with 
opening half-angle (3. The X-ray-emitting gas streams 
along this cone away from the line of centres at speed 
vq (see Fig. Assuming the material is distributed 
symmeUically around the shock cone, the line centroi d 
shift (v) and width are given by (IPollock et all 12004) 



v = — t>ocos/3cos7 (1) 

FWHM = v sin (3 sin 7 (2) 

The angle 7 can be calculated from the orb it 
( Schmutz et all Il997l : Ide Marco & SchmutzL Il999h : 
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Figure 2: The shifts and widths of the emission lines 
in the Chandra-HETGS spectrum of 7 2 Vel, plotted as 
a function of lab wavelength. 



7 = 44° at the time of the Chandra observation. By 
eliminating vq from the above equations, and using 
our measured shifts and widths (v = — 52kms~ 1 , 
FWHM = 1240 km s" 1 ) we find that (3 = 88°. 
Note, however, that the mean observed line shift is 
comparable to the absolute wavelength accuracy of the 
HETGS (~100 km s -1 ; Chandra Proposers' Obser- 
vatory Guide). This means that the measured value 
of v (and hence the derived value of f3) may not be 
completely trustworthy. However, we can say that 
\v\ < 100 km s _1 (since we would unambiguously be 
able to detect a larger value of v). This implies that 
f3 > 85°. 

To understand this value in terms of the dynam- 
ics of the system, we have carried out a hydrody- 
namical simulation of the wind-wind collision. For 
simplicity we have assumed non-accelerating, spher- 
ically symmetric winds. Since there is evidence that 
the O s tar wind collides below its termi nal veloc- 
ity (e.g. ISt-Louis. Willis. & Stevens! Il993h. we have 
adopted an O star wind speed of 1500 km s -1 (cf. 



1990; 



Prinja. Barlow. & Howarthi 
St-Louis et ah. 1993^ We assume the WR 



2300-2400 km s 



star collides at its terminal ve locity (1500 km s 1 ; 
Barlow. Roche. & Aitkeri , 1988b . We adopt mass-loss 
rates of M WR = 1 x 10" 5 M Q yr" 1 and M = 
5 x 10 -7 M Q yr _1 . Figure |4] shows a density map 
from the resulting simulation. 

We have used this simulation to model the X-ray 
li ne emission from 7 2 Vel using the model described 



in iHenlev et alJ d2003h . which calculates X-ray emis- 
sion line profiles from the hydrodynamical simula- 
tion results assuming collisional ionization equilib- 
rium. The resulting profiles are typically blueshifted 
by ~300 km s _1 , in contrast to the essentially un- 
shifted lines in the observed spectrum. 

The discrepancy between the observed and calculated 
line shifts is most likely due to the shock opening half- 
angle in the hydrodynamical simulation being ~40°, 
whereas the simple geometrical model described above 
implies a shock opening half-angle of >85°. In a CWB 
with non-accelerating winds, j3 is a function only of the 
wind momentum ratio (Mwruwr/Mouo)- As this 
ratio increases, [3 decreases, because the more pow- 
erful wind of the WR star overwhelms that of the O 
star. The shock opening angle in our hydrodynami- 
cal simulation therefore depends on the adopted val- 
ues of the wind parameters. However, the shock open- 
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Figure 3: A simple geometrical model for the wind- 
wind collision in a CWB. The black circles denote the 
two stars and the red cone (with opening half-angle f3) 
is the wind-wind interaction region (along which X- 
ray-emitting material is streaming at speed vq). 



ing half-angle inferred from the simple geometrical 
model (J3 > 85°) implies approximately equal wind 
momenta, which is not the case for any sensible set of 
(constant velocity) wind parameters. We therefore sug- 
gest that this large implied opening angle may be evi- 
dence of sudden radiative braking dOwocki & Gayley . 



ZX I u 



z (cm) 



19951: iGavlev. Owocki. & CranmerL I1997TT in which 



the wind of the WR star is rapidly decelerated when it 
encounters the radiation field of the O star. This rapid 
deceleration leads to an increase in the shock opening 
angle (the larger the opacity of the WR star wind to 
the O star radiation field, the larger the shock opening 



angle; see Fig. 3 in lGavlev et all 11997). It also alters 



the Mach number of the wind-wind collision, resulting 
in softer X-ray emission than would otherwise be ex- 
pected. However, as well as affecting the X-ray emis- 
sion, radiative braking may also affect the non-thermal 
radio emission. This is because the orbital variability 
of the non-thermal radio emission (due to free-free ab- 
sorption in the stars' unshocked winds) is dependent 
on the shock opening angle. 

Chapman et al. Jl999h analysed ATCA radio obser- 
vations of 7 2 Vel at four wavelengths using a 
thermal + non-thermal emission model. The ther- 
mal free-free emission is from the unshocked stel- 
lar winds, while the non-thermal synchrotron emis- 
sion is from the wind-wind collision and is attenu- 
ated by free-free abso rption in the unshocked winds. 
Chap man et all (U.999) found that over the range ~2- 



Figure 4: Density map from a hydrodynamical sim- 
ulation of 7 2 Vel. The adopted wind parameters 

are M WR = 1 x 1CT 5 M Q yr" 1 , M = 5 x 
1Q- 7 M Q yr -1 , v WR = vq = 1500 km s -1 . 



4 GHz, thermal and non-thermal emission make ap- 
proximately equal contributions to the observed flux. 
However, their observations were made at two differ- 
ent orbital phases, and so their assumptions of constant 
non-thermal spectral index and constant optical depth 
for the non-thermal emission may not be valid. 

Radio observations with better orbital coverage (S. 
Dougherty, priv. comm.) show no evidence of non- 
thermal emission from Vel. This is presumably be- 
cause 7 2 Vel is a relatively close binary, and so the 
non-thermal emission from the wind-wind collision 
is strongly absorbed by the unshocked stellar winds. 
Also, the population of shock accelerated relativis- 
tic electrons will be quickly reduced by efficient in- 
verse Compton cooling. Thus, it may unfortunately 
be the case that the effects of radiative braking on 
the non-thermal radio emission discussed above will 
be very difficult to observe, b ecause radiative br aking 
is favoured in close systems dGavlev et al. . 1997 ) in 
which any non-thermal radio emission from the wind- 
wind collision will be strongly attenuated by the un- 
shocked winds. 



3 WR140 

The eccentric WC7 
2900 ± 10 days, e 
1990) is another key system in our understanding of 



+ 04.5 binary WR140 (P 



0.84 ± 0.04; Williams et al 
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Figure 5: The co-added first order MEG spectrum of WR140. The spectrum has been binned up to 0.02 A. The 
absorption at long wavelengths is due to the ISM. 



CWBs. Compared with 7 2 Vel, WR140 is hotter and 
about two orders of mag nitude more lum inous. ROSAT 
and ASCA observations JCorcoranl.ll99fih. and m ore re- 
cently RXTE observations ( Po llock et all 120041) show 



metric winds and collisional ionization equilibrium. 



The wind parameters used are from Stevens et al 



a gradual rise in X-ray flux from apastron to periastron, 
followed by a rapid drop in flux at around the time of 
periastron passage. The flux variations agree qualita- 
tively with the theoretical predictions of Ste vens et al 



( 1992) - the increase in flux prior to periastron is due to 
an increase in the density of the X-ray-emitting gas as 
the stars approach each other, whereas the subsequent 
sudden drop in flux is due to increased absorption as 
the line of sight moves from the O star wind into the 
denser WR star wind. 

Figure |5] shows the first-order MEG spectrum from 
a 45-ks Chandra-HETGS observation of WR140 ob- 
tained a few weeks be f ore its most recent periastron 
passage ( Poll ock et all 120041) . As for j 2 Vel, line 
shifts and widths were measured by fitting Gaussians 
to individual emission lines. The measured shifts and 
widths are shown in Fig. |6] plotted against ionization 
potential. The lines are systematically blue-shifted by 
~600 km s _1 , and the line widths are correlated with 
ionization potential. 

As for 7 2 Vel, we have attempted to model the X- 
ray line emi s sion f rom WR140 using the model of 
Henl ey et all (120030 . again assuming spherically sym- 



( 1992). Radiative braking is not an issue in WR140, 
as it is a much wider binary than 7 2 Vel [at the times of 
their respective Chandra observations, the stellar sepa- 
rations were 0.9 A.U. (7 2 Vel) and 4.2 A.U. (WR140)]. 
The calculated line shifts and widths are compared 
with the observed values in Fig. |6] As can be seen, 
there is poor agreement between the predictions of the 
model and the observed shifts and widths. Firstly, the 
observed shifts are smaller than expected. This implies 
that the velocity along the shock cone is lower than ex- 
pected. Indeed, the highest velocity observed in the re- 
solved lines (~2000 km s -1 ) is significantly lower than 
the velocities of the unshocked winds (~3000 km s -1 ; 
Pollock et all 120041) . This may mean that some of 



the energy budget is going into the production of rel- 
ativistic particles. Secondly, the correlation between 
the line widths and ionization potential is not repro- 
duced by the model. This may be indicative of non- 
equilibrium ionization, as it implies that higher excita- 
tion lines (e.g. S XVI) originate further from the line of 
centres where the velocity is larger. This is contrary to 
what is expected if collisional ionization equilibrium 
holds, because the temperature of the shocked gas in 
the wind-wind collision decreases monotonically away 
from the line of centres (hence lines from more highly 
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Figure 6: The shifts and widths of the emission lines in 
the Chandra-HETGS spectrum of WR140 compared 
with shifts and wid ths calculated using the model of 
Henl ev et al.l(l2003h . 



ionized species would be expected to originate nearer 
the line of centres). 

O ver the cour se of eight years at the VLA, 



White & Becker ( 1995) obtained a detailed radio light 
curve of WR140. The observed flux is non-thermal, 
and is highly variable at all three wavelengths observed 
(2, 6 and 20 cm). The radio emission peaks before pe- 
riastron passage, at an orbital phase of (f> ~ 0.7-0.8 
(where <fi = corresponds to periastron). The 2 cm 
emission peaks earliest, and is already on the decline 
while the 20 cm emission is still rising. 



White & Becker ( 1995) were unable to explain the 
variation in the radio flux (which they attributed to 
both a varying intrinsic luminosity and a varying op- 
tical depth) with spherically symmetric winds. In- 
stead they proposed a model in which the wind of the 
WR star is disk-like. The shifts and widths of the X- 
ray emission lines from WR140 are also not well de- 
scribed by our current model with spherically symmet- 
ric winds, as discussed above. This may lead one to 
speculate whether or not they can a lso be explained 
by a d isk-like wind model. However, IWhite & Becker 
Jl995h essentially assume the non-thermal source is 
point-like, whereas theoretical modelling of the ra- 
dio emission from wind-wind collisions shows that 
it is important to consi der extended emission regions 
JPou ghertv et al. . 2003 ). Therefore, it may be possi- 
ble to explain WR140's radio lightcurve with spher- 
ically symmetric winds, without having to invoke a 
disk-like wind model. Furthermore, the variation in 
the observed X-ray flux is in good qualitative agree- 
ment with the predicti ons of a spheric a lly symmetric 
collid ing wind model dCorcorari Il99d : IPollock et all 
1999). It seems that more detailed modelling of the 
X-ray emission lines and the radio emission, and more 
spectral information from the radio are required to re- 
solve the issue of the nature of WR140's winds. (See 
also Pittard et al., these proceedings.) 



4 Conclusions 



High-resolution X-ray spectra allow us to probe the 
structure and dynamics of the wind-wind collision in 
CWBs. In the two key systems described here, the X- 
ray spectra have revealed evidence of a variety of in- 
teresting phenomena, such as sudden radiative braking 
and non-equilibrium ionization. From a radio point-of- 
view, the former is of particular interest, as radiative 
braking increases the shock opening angle, which in 
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turn influences the variability of the non-thermal emis- 
sion. 

Future modelling of the X-ray emission lines needs 
to take into account the effects of the stars' radiation 
fields and the acceleration of relativistic particles on 
the hydrodynamics, as well as calculating the ioniza- 
tion balance self consistently (rather than assuming 
collisional ionization equilibrium). Furthermore, while 
line emission is generally considered a thermal pro- 
cess, ions may equally well be excited by collisions 
with non-thermal electrons. If the relative abundance 
of non-thermal electrons is large (the energy density of 
relativistic electrons in CWBs may b e as large as ~1% 



of th e thermal energy density; see Doug herty et al 



2003) these too will have to be included in the mod- 
els. 
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